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Abstract 

We study the decays of Bg d l~^l~'y (/ = e, /i,r) within the hght-front modeL We 
calculate the tensor type form factors and use these form factors to evaluate the decay 
branching ratios. We find that, in the standard model, the branching ratios of -Bs(d) 
l+l-j (/ = e,/i,r) are 7.1 x lO'^ (1.5 x 10"^°), 8.3 x 10"^ (1.8 x 10-^°), 1.6 x 10"^ 
(6.2 X 10~^°), respectively. 



1 Introduction 



It is well known that B physics is important to determine the elements of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix[0 and physics beyond the standard model. Recently, 
the interest has been focused on the rare B meson decays induced by the flavor changing 
neutral current (FCNC) due to the CLEO measurement of the radiative 6 — 57 decay 0] . 
In the standard model, these rare decays occur at loop level and provide us information 
on the parameters of the CKM matrix elements as well as various hadronic form factors, 
such as the -B meson decay constant fs- 

As in the decays of B^ l^ui , the helicity suppression effect is also expected in 
the flavor changing neutral current processes of -Bs,d l^l^ . These decays are sensitive 
probes of top quark couplings such as the CKM elements Vts{d)- The decay widths of 
these leptonic decay modes are given by: 

r(s, - i^n = ^^gpi < (£1) \v,,v:,\-cl , (1) 

where Gi? is the fermi constant, Mb^ and mi are Bq meson and lepton masses and Cio 
is the Wilson coefficient. Form Eq.(|l|), one has that B{Bs e~^e~ , fi~^ fi~ ,t^t^) ~ 
(6 X 10-^2.6 X 10-\1.0) X 10~^ and B{Bd e+e- , fi+ fi' ,t+t-) ~ (4.2 x 10"^ 1.8 x 
10-2,4.5) X 10-s by taking \VtbV*\ = 0.04, iVtbV;^] = 0.01, /b, ^ ~ 200 MeV, tb, ~ 1.61 ps 
and tbj — 1.5 It is clear that the rates for the light lepton modes are too small to 

be measured due to the helicity suppressions, while that for the r channel, although there 
is no suppression, it is hard to be observed experimentally because of the low efficiency. 

It has been pointed out that the radiative leptonic B decays of B^ l^i^n = 
e, fi) have larger decay rates than that purely leptonic ones [^, |, ^ ^ |10|, |TT|, . 
Similar enhancements were also found for the radiative decays of -B^.d l^l^'y in the 
constituent and light-cone sum rule quark models |1^, |T^. In fact, the decay amplitudes 
of Bs^d l^l^'~f can be divided into the "internal-bremsstrahlung" (IB) parts, where 
the photon emits from the external charged leptons, which are still helicity suppressed 
for the light charged lepton modes, and the "structure-dependent" (SD) ones, in which 
one of the photon comes from intermediate states of Bq — > which are free of the 

helicity suppression. Therefore, the decay rates of Bq l^l^'-f { I = e, fi ) might have 
an enhancement with respect to the purely leptonic modes of Bq if the SD 

contributions to the decays are dominant. 
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In this paper, we will use the light front quark model [1^, |TB|, |TB[ to evaluate 

the hadronic matrix elements and study the decay rates for i?^,^ l^l~'~f. It is known 
that as the recoil momentum increases, we have to start considering relativistic effects 



seriously. The light front quark model [|T^] is the widely accepted relativistic quark model 
in which a consistent and relativistic treatment of quark spins and the center-of-mass 
motion can be carried out. In this work, we calculate for the first time the tensor form 
factors in P — > 7 (P is a pseudoscalar meson) directly at time-like momentum transfers by 
using the relativistic light-front hadronic wave function. Within the framework of light- 
front quark model, one can calculate in the frame where the momentum transfer is purely 
longitudinal, i.e, p± = and p"^ = p'^p^ , which covers the entire range. Thus, we give their 
dependence on the momentum transfer in whole kinematic region of < < p^ax- 

The paper is organized as follows. In Sec. 2, we present the decay amplitudes of 
Bs,d ~^ l^l~'~f. We study the form factors in the P ^ 7 transition within the light front 
framework in Sec. 3. In Sec. 4, we calculate the decay branching ratios. We also compare 
our results with those in literature [0, |T9|. We give our conclusions in Sec. 5. 

2 Decay Amplitudes 



The main contributions for the processes of Bq l^l 7 (/ = e, yU, r) arise from 



the effective Hamiltonian that induces the purely leptonic modes of Bq 
QCD-corrected amplitude for h l^l^q in the SM is given by I^U], pT| : 



/+/-. The 



M 



2—qia^up''{mbPii + mqPL)bl-f^l 



(2) 



where q = d or s, Pl,r = (1 T 75)/2, and p is the momentum transfer, which is equal 
to the momentum of the lepton pair. The Wilson coefficients of Cj, Cg^^ and Ciq can be 
found, for example, in |^ ^ respectively. 
The amplitude for Bq can be written as: 

M{Bq ^ /+/-7) = MiB + MsD (3) 

where Ai ib and A4sd represent the IB and SD contributions, respectively. For the IB 
part, the amplitude is clearly proportional to the lepton mass m/ and it is found to be 



M 



IB 



GpQ. 

-ie^=-VtbVt* fs, Ciomi 
V 27r 



2pi ■ q^ 2p2 ■ q^ 



(4) 
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where Pb, Pi, P2 and q^, are the momentum of Bq, /+, / and 7, respectively. In Eq.(4), 
the Jb^ is the Bg meson decay constant, defined by: 



(5) 



When a photon emitted from the charged internal line, the amplitude is suppressed by 
a factor mf/m^ in the Wilson coefficients and the main contribution comes from the 
photon radiating from the initial quark line. Therefore A4sd can be written as 



M 



SD 



+ Cio{^{qj)\q^^,PLb\B{p + q^))h^^^l 



-{l{q^)\qi(J^r.p''PRh\B{p + q^))lYl\. 



(6) 



From the amplitude in Eq.(6), we see that to find the decay rates, one has to evaluate 
the hadronic matrix elements. These elements can be parameterized as follows: 

FvA 



{l{<li)\(lli^lbh\Bg{p + q^)) = 
{l{(li)\Q.lABqip + (li)) = 

{-l{q^)\qi(r^uP''lbh\Bq{p + q^)) 
{-i{q^)\qi(J^yp''h\Bq{p + q^)) 



Mb, 



W *a 3 u 



ieFr 



TA 



{P-Qy 



eFTv£nuai3e*'"q^P^, 

where the is the photon polarization vector, q^ and p + q-y are the four momenta of 
the photon and the Bg meson, and Fva, Fyv, Ft a and Ftv are the form factors of axial- 
vector, vector, axial-tensor and tensor, respectively. Form Eqs.(7) and (8), we rewrite the 
amplitude of Eq. (6) as. 



(7) 



(8) 



M 



SD — 



2V27r 



VM e^.a/3e* A 1^1 + Cin.l 



+ i 



{p ■ q,)el - {e* ■ J [b 1^1 + D l^^.l 



(9) 



where the factors of A-D are defined by 



A = 


r^ff 


-2C,^FrA(p'), 

pi 




B = 


% ^-(-') 


-2Cr^FMp'), 




C = 








D = 






(10) 



respectively. The form factors of Fva and Fyv have been evaluated previously in the 
light-front model [|^, while that of Ft a and Ftv shall be studied in the next section. 



3 Form Factors in the Light Front Model 



In this section, we will use the light-front approach to calculate the tensor type form 
factors in Eq. (|^) for Bg 'y {q = s or d) transition. In this approach, the B meson 
bound state consists of an anti-quark b and a quark q with the total momentum of [p + q-y) 
and it is given by 

\B{p + q^)> = j [dh][dk2]2(2T^f5''{p + q^-h- k^) 



A1A2 



X k^)hUki. XM{k2, \2)\0 > , (11) 



where fci(2) is the on- mass shell light front momentum of the internal quark b{q), the light 
front relative momentum variables {x, k±) are defined by 

kf = x{p + , ki± = x{p + g^)_L + k^ , (12) 

and 

/ 2k^k~^ \ ^ 

^'B'''{x,k^)=[ ' ' ^ u{k,,X,)j'v{k2,X2)4>{x,k^), (13) 

\Mq — [mq — mb) ) 

with 0(a;, fcj,) being the momentum distribution amplitude. The amplitude of 0(x, A;_l) 
can be solved in principle by the light-front QCD bound state equation P^, Here, we 
use the Gaussian type wave function: 



(.,.,)^iv/§exp -±,1. (14) 



where 



. „ , dk+dk_i_ , / TT V 



2/ 2Mo " X 1-x ' 

Y: u{k, XMk, A) = ^ , E XMk, X) = . (15) 
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We note that the wave function in Eq.(14) could be also applied to many other of hadronic 
transitions. For the gauged photon state, one has 



|7(g^) > = iV'|a+(g^, A) + ^ / [dh][dk2]2{2T:f5\q, -k,- k^) 
I A1A2 



X ^]f'\q,. ku k2)K{ku Ai)4(fci, A2) Mo > , (16) 



where 



cr_L ■ kx, +imi^ \ 

— o-± ■ e± > Xxi , (17) 



with 



ED 



^7! + ^2x + ^ 



q+ kt kt ^ ' 

In Eq. (|T^, we have used the two- component form of the light-front quark fields |T^ . 



Since the physically accessible kinematic region is < < p^^^^ where p^nax = 
to calculate the matrix elements in Eq. (^, we choose a frame where the transverse 
momentum p_L = 0. Then = p^p^ > which can cover the entire range of the 
momentum transfers. By considering the "good" component of /i = -|-, the tensor current 
in Eq. (|^) can be rewritten as: 

< 7(g7)l(9i7°75&- - gi7%5&+)|5(p + g^) > = cFta (el ■ g^x) > 

<liq,)\iqtl'b_-q^^%^)\Bip + q,)> = -teFrve^'^qy, , (19) 

where g+(6+) and q^{h_) are the light-front up and down component of the quark fields. 
In the two-component form [1^ , they are expressed by 



= ( S ) (20) 

and 

= . a, + /3,„,),, = ( _^ ) . (21) 

In Eq. (21), Xq is a two-component spinor field and cr is the Pauli matrix. The form 
factors of Ft a and Ftv in Eqs. (|T^ are then found to be: 

dx d'^k I 



ji A^ + A^kie i B, + B,kie \ 



2 (27r)' 

[ 1 Ci + C2 fcie 1/^1 + 1^2 A;i0 ' 

X < s r^F h -- 



3 ml + P, 



3 171^ + k'j_ 



—p^ ^— -| (x' + x - 2x'x) [x{x - 1) - a;(2a; - 1)] ki 

xx''^[i — x')[i — x) y 

+x [{x — x') + 2x'x{l — x)] ml + 2x^(1 — x')^mbmq^ , 

— Aix' + x- 2x'x)(l - 2x)kl 

xx'^{l-x'){l-x)y 

+2x{l — x')mqmb — x{l — 2a;) (1 — a;')^m^ + x'{l + x'x — 2xx''^)m 
2 



x'x{l — x){l 



— \ {x' + x- 2x'x){l -2x + 2x^ - x'x)kl 

— X ) 



+2xx'{l - x'){l - x)mqmb + (1 - x'){x' + x- Ax'x + 2xV)mJ 

2{x- x') 



—^l-(x' + x-2x'x)(l-2x)kl 

x'x{i-x){i-x'y\ ^ j\ J ± 

— (1 — 2x)x''^{l — x)ml + (1 — x'){x' + X — Sx'x — 2x^(1 — x))m 
xx''^{l-x'){l-x) 



{x — x' -\- xx'){x' + X — 2x'x)k^ 



-2x^(1 - x'^mbTUq + x' Ux - x') - 2(1 - x')x^ 
2{x- x') 



xx'\l-x'){l-x) 



x' + X — 2x'x)kj_ + x'{l — 2x + xx')mf^ 



—mgX{l — x') — 2x{l — x')mq > , 



(1 -x){l-2x + x'){x' + x- 2x'x)k'^ 



x{l — x)x'{l — x'Y 



(1 - x'){x' + x-2x^ - 2x'x + 2x'^x')ml + 2x''^(l - xfrnqTUb 



2{x- x') 



x'{l — x)x{l — x') 
+(1 — x'){x' — X — x'x)m'i > , 



— Ux' + x- 2x'x)k'i + - x)ml 



N 



2x(l - x) 
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1/2 



- {ruq - mbf 



I dkz 
dx 



exp 



2a;2 I ' 
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^ 1 d<l>{x,kl) 
dkl ' 

x' = x(^l-^y k = {k^X)- (24) 

To illustrate numerical results, we input uiq = rris = 0.4, rrib = 4.5 in GeV, and 
u = 0.55 GeV which is determined by fitting Jb, = 200 GeV via Eq. (5). The values of 
Ft A and F^v in the entire range of p'^ are shown in Fig. 2. We note that the tails of Ftv,ta 
at the large momentum transfers go down may be the long distance contribution associated 
with B — B* — 'J vertex, which is not included in the present work. It is interesting to 
note that the formulas in Eqs. (0) and (^) can be used for other pseudoscalars, such 
as pions and kaons, to the photon transitions as well once we put in the corresponding 
masses. 



4 Decay Branching Ratios 



The partial decay width for Bg l^l 7 in the B rest frame is given by 

dV = -J—\M\^(2n)^5\PB - Pi - P2 - gJ , "jf , "(f' , "^f' , (25) 
2Mb \ b fi f2 '^T^(27r)32E^ (27r)32^i (271)32^2 ^ ' 

where the square of the matrix element can be written as 

I-MT = \MiB? + \MsD? + 2Re{MiBM*so) ■ (26) 

To describe the kinematic of Bg we defined two variables of x-y = 2Pb ■ q^/Mb 

and y = 2Pb ■ Pi/Mb- One can easily write the transfer momentum p'^ in term of x^ as 

/ = M|(l-x,). (27) 

The double differential decay rate is found to be 



dXydX 2567r3 



^-|M|' = Cp(x„A), (28) 



where 



C = al^^l'GlMl, , (29) 



and 



p{x^,\) = pib{x^,X) + Psd{x^,>^) + Pin{x^,\), (30) 
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with 



PiB 



PSD 



Pin 



2 



' xt^- 2x^ + 2- Ari 



2n(— + 



A2 (1-A)2' 



M 



x'l {\A\' + \B\') (1 - + 2n) - 2(1 - x^){\ - y) 



+ (|C|2 + (1 - - 2ri) - 2(1 - a;^)(A - 



+2Re{B*C + A*D){1 - x^)(2A - 1) 



+ Re{D) 



x^{l - 2A) 



(31) 



A(l-A) ' ' A(l-A) j ■ 

Here A = {x-y + y — l)/x-y and n = mf/M'^ and the physical regions for x^ and A are given 
by: 



< a;^ < 1 - 4r/ . 



2 2V 



4n 



1 — x^, 



< A<^ + ^Wl 



4r/ 



1 — 



(32) 



For the Wilson coefficients Cj and Cio, we use the results given by Refs. [|20|, |21[] and 
they are 

C7 = -0.315 , Cio = -4.642 . 

The analytic expressions for Cg^'^ in the next-to-leading order approximation is given 
bylH: 



with 



Cg + 0.124w;(s) + g{m„ s)(3Ci + C2 + 3C3 + C4 + 3C5 + Cg) 
-i(7(^g, s){C3 + 3C4) - i^/(m6, s){4Cs + 4C4 + 3C5 + Cg) 

+ ^(3^3 + ^4 + 3^5 + ^6), 



(33) 



C9 

(3Ci + + 3C3 + C4 + 3C5 + Cq) 

(C3 + 3C4) 
(4C3 + 4C4 + 3C5 + C6) 
(3C3 + C4 + 3C5 + Ce) 



4.227 , 
0.359 , 
-0.0659 , 
-0.0675 , 
-0.00157 



(34) 
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and s = p'^/ml. In the Eq.(^), the function of w{s) comes from the single-gluon correction 
to the matrix element of Og and its form can be found in Refs. , while that oig{rhi, s) 

from the one loop contributions of the four quark operators Oi - Oe, given by pll, p3|: 



g{mi,s) 



-- Inmj H 1 — Vi 

9 27 9^^ 
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X 



(lnl±^ 
2 arctan 



1 



+ 

for yi<l 
for yi> I 



(35) 



with yi = rhi' / s and rhi = rrii/mb. By taking into account the long distance effects mainly 
due to the J/ip family resonances, one may use the replacement] 

rnyBr{V 



712 2^ 



(36) 



where rriy = mv/rrih and ttot = T/^^fe • The masses and decay widths of the corresponding 
mesons in Eq.(35) are listedQ] in table 1. 



Table 1: Charmonium mass and width 



Meson 


Mass (GeV) 


Br{V l+l~) 


Vtotai{MeV) 




3.097 


6.0 X 10~2 


0.088 




3.686 


8.3 X 10-3 


0.277 



In Figs. 2 and 3 we present the differential decay rates of Eg fi~^fi~j and Bg t^t~j 
as functions of x-y, with and without resonance {J/ip (^f^d ip') contributions. From these 
figures we see that the contributions from the IB parts, corresponding to small x-y region, 
are infrared divergence. To obtain the decay width of Bs t^t^^, a cut on the photon 
energy is needed. Our results for the integrated branching ratios without and with long- 



distance effects as well as those given by the constituent quark model |14] and the light- 



cone QCD sum rule model |T^, [T^ are summarized in tables 2 and 3, respectively. Here, 



we have used the cut value of 5 = 0.01 and rrir = 1.5 GeV, m, = 0.4 GeV, W-fhV^ 



tbVtsl 



0.045, 



tb^tdl 



0.01, t{Bs) = 1.61 X 10-12 5 and T{Bd) = 1.5 x lO^^^ ^||_ 
We now compare our results with those in the literature [|T3], 0, As shown in 
Table 2, the decay branching ratios for Bg l^l^^ (/ = e, /i) without long distance 
contributions found in our calculations are smaller than those in the constituent quark 
model H^, whereas for B^i —>■ the statement are much smaller than ones. It is 

mainly due to that in the constituent quark model the decay rate of Bg 
is proportional to the inverse of the quark mass m^. It is interesting to see that our 
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Table 2: Integrated branching ratios for the radiative leptonic Bs^d decays without long- 
distance 



Integrated Decay 


IB 


SD 


IN 


S^nm 

kJ Lllll 


Ref. 

fl4l 
[I'll 


Ref. 

13 IQl 


10^^(5, ^ e+e-7) 


3.1 X 10-^ 


3.2 


5.2 X 10-*^ 


3.2 


6.2 


2.35 


1095(5, ^ /i+/i-7) 


5.5 X 10-1 


3.2 


8.2 X 10-2 


3.8 


4.6 


1.9 


10^B{Bs T+T--i) 
(Cut 6 = 0.01) 


13.0 


0.6 


6.4 X 10-1 


14.2 




9.54 


W^°B{Bd e+e-7) 


1.4 X 10-^ 


1.5 


2.4 X 10-6 


1.5 


8.2 


1.5 


10'^B{Bd ^ //+/^-7) 


2.6 X 10-1 


1.5 


3.8 X 10-2 


1.8 


6.2 


1.2 


10i°5(5d ^ r+r-7) 
(Cut 5 = 0.01) 


5.7 


2.4 X 10-1 


2.7 X 10-1 


6.1 







Table 3: Integrated branching ratios for the radiative leptonic B^ d decays with long- 
distance 



Integrated Decay 
Branching Ratios 


IB 


SD 


IN 


Sum 


Ref. 

|19 


10^5(5, ^ e+e-7) 


3.1 X 10-^ 


7.1 


5.8 X 10-6 


7.1 




10^5(5, ^ /i+/i-7) 


5.5 X 10-1 


7.1 


9.4 X 10-2 


8.3 




10^B{Bs r+r-7) 
(Cut 6 = 0.01) 


13.0 


1.9 


7.9 X 10-1 


15.7 


15.2 


10i°5(5d ^ e+e-7) 


1.4 X 10-5 


1.5 


2.4 X 10-6 


1.5 




10105(5, ^ /i+/i-7) 


2.6 X 10-1 


1.5 


3.8 X 10-2 


1.8 




10i05(5rf r+r-7) 
(Cut S = 0.01) 


5.7 


2.3 X 10-1 


2.5 X 10-1 


6.2 





results are close to those in the light-cone QCD sum rule model |[13| , p!9| . We note that 
the SD contributions for the decays in both constituent quark and light-cone QCD sum 
rule models are sensitive to the values of the decay constants fs^ ^• 

5 Conclusions 



We have studied the decays of 5, ^ l^l~'y within the light-front model. We have 
calculated the tensor type form factors and used these form factors to evaluate the de- 
cay branching ratios. We have found that, in the standard model, the branching ratios of 
Bs(d) e+e-7, Bs(d) /i+yu"7 and Bs(d) ^ r+r-7 are 7.1 x 10-^ (1.5 x 10-i°), 8.3 x 10"^ 
(1.8 X lO-i*') and 1.6 x 10-^ (6.2 x lO-i'^), respectively. Comparing with the purely lep- 
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tonic decays of Bg^d l^l^, we have shown that the branching ratios of Bg^d l^l^'j 
have the same order of magnitude for the /i channel but that of B^ d ~^ e"'"e~7 are much 
larger. We conclude that some of the radiative leptonic decays of Bg^d could 
be measured in the B factories as well as LHC-B experiments, where approximately, 
6 X 10^^(2 X 10^^) Bd{Bs) mesons are expected to be produced per year. 
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Figure Captions 



Figure 1: The values of the form factors Ft a (sohd curve) and F^v (dashed curve) as 
functions of the momentum transfer for Bg 7. 

Figure 2: The differential decay branching ratio dB{Bs — > jJL^ jjT ^) / dx^ as a function of 
x-y = 2E^/Mb^ with(solid curve) and without (dashed curve) long distance. 

Figure 3: The differential decay branching ratio dB{Bs T~^T~j)/dXj as a function of 
x-y — 2Ej/Mb^ with (sohd curve) and without (dashed curve) long distance. 
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Figure 1: The values of the form factors Ft a (sohd curve) and Ftv (dashed curve) as 
functions of the momentum transfer for — > 7. 
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Figure 2: The differential decay branching ratio dB^Bg j)/dxj as a function of 

Xj — 2Ej/Mb^ with (sohd curve) and without (dashed curve) long distance. 
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